PITHA 08/01 



Determining the CP parity of Higgs bosons 

at the LHC 
in their t decay channels 

Stefan BergeH, Werner BernreutheiHI, Jorg Ziethe 

Institut fiir Theoretische Physik, RWTH Aachen, 52056 Aachen, Germany 



Abstract 

If neutral Higgs bosons will be discovered at the CERN Large Hadron Collider (LHC) then 
an important subsequent issue will be the investigation of their CP nature. Higgs boson 
decays into x lepton pairs are particularly suited in this respect. Analyzing the three charged 
pion decay modes of the x leptons and taking expected measurement uncertainties at the 
LHC into account, we show that the CP properties of a Higgs boson can be pinned down 
with appropriately chosen observables, provided that sufficiently large event numbers will 
eventually be available. 
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Introduction: The major physics goal of the upcoming experiments at the CERN Large 
Hadron Collider (LHC) is the exploration of the hitherto unknown mechanism of electroweak 
gauge symmetry breaking which, in the context of the standard model of particle physics 
(SM) and many of its extensions, is tantamount to searching for Higgs bosons, spin-zero 
and electrically neutral resonances with masses of (a few) hundred GeV (see, e.g., [1,2] for 
reviews). If (one or several types of) Higgs bosons are found then the next issue will be 
the determination of their properties - in particular their parity (P) and charge conjugation 
times parity (CP) quantum numbers, respectively, which yield important information about 
the dynamics of these particles. While the SM Higgs boson is parity-even, SM extensions 
also predict parity-odd state(s) or, if the Higgs boson dynamics violates CP, states of un- 
defined CP parity with Yukawa couplings both to scalar and pseudoscalar quark and lepton 
currents. Higgs sector CP violation is a fascinating possibility, especially in view of its po- 
tentially enormous impact on an important issue of the physics of the early universe, namely 
baryogenesis [3]. These new interactions can be searched for at the upcoming generation 
of colliders in several ways (see [2] for a collection of recent results). The decays of Higgs 
bosons to x^x+ leptons and/or - if the Higgs bosons are heavy enough - to tt quarks are par- 
ticularly suited in this respect [4-10]. In this letter we show that, at the LHC, the CP nature 
of a neutral Higgs boson can be pinned down with appropriately chosen observables in its 
x^x+ decay channel where the x decay into three charged pions, provided that sufficiently 
large event numbers will eventually be available. 

The analysis of this letter applies to any neutral Higgs boson hj with flavor-diagonal cou- 
plings to quarks and leptons / (with mass nif) 

Ly = ~{V2Gfy/^Y.'^fiajfff + bjffiy5f)hj , (1) 

j,f 

where Gp is the Fermi constant and ajf and bjf are the reduced scalar and pseudoscalar 
Yukawa couplings respectively, which depend on the parameters of the (effective) Higgs po- 
tential of the respective model. In the SM, af = I and bf = 0. As far as SM extensions 
are concerned we consider here, for definiteness, models with two Higgs doublets, such 
as the non-supersymmetric type II models and the minimal supersymmetric SM extension 
(MSSM) (see, e.g., [1, 2]). These models contain three physical neutral Higgs fields hj in 
the mass basis. If Higgs sector CP violation (CPV) is negligibly small then the fields hj 
describe two scalar states h,H (bf = 0) and a pseudoscalar A {af = 0). In the case of Higgs 
sector CPV, the hj have non-zero couplingil] ajf and bjf to quarks and leptons which lead 
to CP-violating effects in hj — > // already at the Born level [4]. This is in contrast to the 
couplings to W+VK^ and to ZZ boson pairs of such a state of undefined CP parity. At the 

'They can be parameterized in terms of the ratio of the Higgs field vacuum expectation values tanp = V2/V1 
and a 3 X 3 orthogonal matrix that describes the mixing of the neutral spin-zero CP eigenstates [11]. 
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Born level, only the CP = +1 component of hj couples to W*"!^" and to ZZ. The coupling 
of the pseudoscalar component of hj - if there is any - to VF+V7^ and to ZZ is likely to be 
very small as it must be induced by quantum fluctuations. Thus, the observation of Higgs 
boson production in weak vector boson fusion W^W ,ZZ — » hj or of the decay channels 
hj W^W ,ZZ would tell us that hj has a significant scalar component. However, the 
question would remain whether or not hj is a pure J^*^ = 0++ state or if it has a significant 
pseudoscalar component, toc§ . This can be answered by investigating the x decay channel 
of this particle. 

X spin observables: In the following, we choose the generic notation (|) for any of the 
neutral Higgs bosons hj discussed above. The observables discussed here for determining 
the CP quantum number of (|) in its decay channel (|) x^x+ may be applied to any Higgs 
production process. At the LHC, this includes the gluon and gauge boson fusion processes 
gg ^ ^ and qjqj ^q'jq'j, respectively, and the associated production tt^ or bb^ of a light 
Higgs boson. We consider the following semi-inclusive reactions 

/ ^ ^ x-(kx,a) +x+(kx,(3) +X, (2) 

where / is some partonic initial state, and kf = — k^ are the 3-momenta of x^ and x+ in 
the XX zero-momentum frame (ZMF), and a, |3 are spin labels. Here we make use of the fact 
that at colliders polarization and spin correlation effects are both measurable and reliably 
predictable for tau leptons. 

Let us assume that experiments at the LHC will discover a neutral boson resonance in a 
reaction of the type ^ and a sufficiently large sample will eventually be accumulated. The 
spin of (|) may be inferred from the polar angle distribution of the tau leptons. Suppose 
the outcome of this analysis is that (|) is a spin-zero (Higgs) particle. One would next like 
to determine its Yukawa coupling(s), and specifically like to know whether (|) is a scalar, 
a pseudoscalar, or a state with undefined CP parity. This can be investigated by using the 
following CP-even and -odd observables involving the spins of x, x, and we emphasize that 
all of them should be used. The correlation resulting from projecting the spin of x onto the 
spin of X, 

^=St-Sx, (3) 

is the best choice for discriminating between a CP = ±1 state [5]. Here StjSf denote the 
spin operators of x^ and x+, respectively. This can be understood as follows. If (|) is a scalar 
(7^^ = 0++) then xx is in a ^Pq state, and (sx • Sx) = 1 /4. If (|) is a pseudoscalar (7^^ = 0^+) 



^ Nevertheless, it can be checked with appropriate correlations whether or not there is a sizeable effective 
pseudoscalar coupling to weak vector bosons [12]. 
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then XX is in a ^5o state and (Sx ■ Sf) = — 3/4, which is strikingly different from the scalar case. 
For general couplings ([B one gets (5) = (a^ — 3&x) / +4^?^), using that m^^ lUx [5]. 
If y^p = —a^bx 7^ 0, the Yukawa interactions of (|) to x leptons are not CP-invariant. This 
leads to CP- violating effects in the reactions Q. For an unpolarized initial state i, a general 
kinematic analysis of Q yields the following result [4, 6] . If C-violating interactions do not 
matter in ^ then Ly (which is C-invariant, but P- and CP-violating) induces two types of 
CPV effects: a CP-odd x^x+ spin- spin correlation and a CP-odd x polarization asymmetry 
which correspond to the observables 

Sep = kx ■ (sx X Sx) , S'cp = kx ■ (sx - Sx) • (4) 

Here, kx = kx/|kx| in the XX ZMli. The CP-odd and T-odcf variable Sep measures a cor- 
relation of the spins of the x and x+ transverse to their directions of flight. A non-zero 
expectation value is generated already at tree level, {Sep) = —a-cbz/ {a^ + b\) [4], which can 
be as large as 0.5 in magnitude. The variable measures an asymmetry in the longitudinal 
polarization of the x^ and x+. As it is CP-odd but T-even, a non-zero (S^p) requires both 
icp 7^ ^ ^'^^ ^ non-zero absorptive part of the respective scattering amplitude. 
(Multi) pion final states: The polarization and spin-correlation effects induced in the 
x^x+ sample from Higgs boson decay lead in turn, through the parity-violating weak decays 
of the X leptons, to specific angular distributions and correlations in the respective final state 
which can be measured with appropriately constructed observables (see below). In order to 
obtain sufficient sensitivity to the CP properties of the Higgs boson resonance, one should 
consider x decay channels that have both good to maximal x-spin analyzing power and allow 
for the reconstruction of the x decay vertex, i.e., the x rest frame, which is essential for an 
efficient helicity analysis. 

We recall the x-spin analyzing power of the final state a in the decay x^ ^ a + Vx, that is, the 
coefficient Ca in the distribution F^^ JFa/JcosGa = (1 -\-Ca cos0a)/2, where 0a is the angle 
between the x^ spin vector and the direction of a in the x^ rest frame (c.f., e.g., [13]). CP 
invariance which is, as known from experiments, a good symmetry in x decays at the level 
of precision required here implies that the x-spin analyzing power of a in x+ — a + Vx is 

The 71^ channel is known to have maximal spin-analyzing power, c^- = 1 . In the case of, 
e.g., multi-pion decays of the x, this optimum analyzing power can also be achieved if all 
the pions are observed and the dependence of the hadronic current on the pion momenta is 



^Two more terms can appear in the squared matrix element of They are obtained by replacing kx — > p 
in (|4|, where p is the direction of one of the colliding beams in /. However, for resonant (|) production only the 
observables (|4]i are of interest. 

'^Here T-even/odd refers to a naive T transformation, i.e., reversal of momenta and spins only. 
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known [14-16]. The latter is obtained using empirically tested matrix elements and fits to 
measured distributions. 

In the following section we consider the decay n^n^n'^Vi and the corresponding decay 
of x+. As this decay proceeds to a large extent via the ai resonance, we use for 
the description of the three-pion final state. The measured pion momenta in the laboratory 
frame allow, using known kinematic distributions [13-15], the separation of the longitudinal 
{an) and transverse (air) helicity states of the cii. This leads to an optimal spin analyzing 
power, c = ±1 for an and air, respectively. Moreover, the measured pion momenta yield 
the X decay vertex and, in turn, the x rest frame (see below). 
We thus investigate reactions of the following type: 

^x-(k„a)+x+(kx,(3)+X ^a(qi)+Mq2)+X, (5) 

where a and b denote here the a^ and a^ resonances, respectively. Without loss of generality, 
we consider the and a^^ states and comment on how the resulting distributions change 
when one or two air are involved. The momenta of x~ and x"^ in ([5]) are defined as above in 
the XX ZMF, while the momenta qi and q2 refer to the x^ and x+ rest frames, respectively. 
Let's now come to the equivalents of the 5,- at the level of the final states a,b. The spin 
correlation (5) leads to a non-isotropic distribution in coscp, where cp = Z(qi,q2). If no 
phase space cuts are applied - modulo cuts on the invariant mass Mxf of the x pair - this 
opening angle distribution is of the form [5]: 

1 d^^ab - _z)^- coscp), D^-, = ^CaCb{s,-S,). (6) 



o^^Jcoscp 2' ""3 

The coefficients D- + = + are 0.33 and —1 for the channels (b (0++) — * XX — ayai ,■ 
and (|)(0^+) ^ XX ^ aua\j, respectively, if ij = LL,TT, while they change sign if ij = 
LT , TL. Thus, for ij = LL.TT the ai momenta qi, q2 are predominantly parallel in the case 
of a pseudoscalar (|), while for a scalar ^ they tend to be antiparallel. 

The equivalents of the CP-odd spin observables Sep and 5^^ at the level of the final states 
© are [5,6]: 

Ocp = (kT-kx)-(q2xqi)/2, O^p = kx ■ qi -kx ■ % ■ (V) 

In general, ifb^a the averages of dV]) should be taken for events ^ plus the charge conju- 
gated events db. Asymmetries corresponding to ([7]) are: 

KdO>0)-N.dO<0)^ (8) 

where A^^^ is the number of events in the reaction ([5]). If no phase-space cuts, besides cuts on 
Mxx, are imposed then [6] 

4 2 

{Ocp)ab = -g(^aC-b {Scp) , {Ocp)aa = {S'cp) , 
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A{Ocp) = TT{Ocp)ab^ Moip) = {Ocp)aa- 



(9) 



Here the relations involving O^p are given for simplicity for the diagonal channels ad only. 
The observable Ocp measures the distribution of the signed normal vector of the plane 
spanned by qi,q2 with respect to the direction of flight. If y^p ^ then this distribu- 
tion is asymmetric. If (|) were an ideal mixture of a CP-even and -odd state, \ax\ = \hx\, the 
asymmetry corresponding to Oqp would take the value \A{Ocp) \ = 0.4 in the auaij channels 
(/, j = L,T). As already mentioned, (O^p) ^ also requires, besides y^p ^ 0, an absorptive 
part in the scattering amplitude z ^ (|) ^ x^x+, which is small in our analysis below. There- 
fore, we do not consider this observable any further. 

Results: For non-standard Higgs bosons (|) and large tan (3, the associated production with 
bottom quarks, gg bb^, is considered to be the most promising mode in the (|) — XX decay 
channel at the LHC [17, 18]. The results shown below are not only applicable to this produc- 
tion process, but also to gluon-gluon fusion, gg — > (|), and vector boson fusion. The reason is 
that our normalized distributions do not depend on the (|) momentum if no detector cuts are 
applied. Furthermore we show for (|) ^ xx ^ aidi that detector cuts have only a very small 
effect on these distributions for Higgs masses larger than 200 GeV. Thus, our results will not 
change significantly for the different Higgs production modes or if initial-state higher-order 
QCD corrections are taken into account. We have, therefore, computed in this analysis all 
distributions for a generic 2^1 Higgs boson production process at leading order. 
As emphasized above, the determination of the distributions of cos cp and of the observable 
Ocp requires the reconstruction of the x^ rest frames. For the decay channels ©, the of 
momenta in the laboratory frame and the x^ decay vertices can be obtained from the visible 
tracks of the three charged pions. The x^ production, i.e., the Higgs production vertex, can 
be reconstructed from the visible tracks of the charged particles/jets produced in association 
with the (|) [19]. Using that, for each x, = q'^^ + qi^, = — k^, = (the tilde refers 
to the laboratory frame), and k-c = Kx, where x is the unit vector along the line connecting the 
X production and decay vertex, the factor K is obtained by solving this system of equations. 
For each x lepton, we obtain two solutions, and we select the solution for which the sum of 
the transverse x^ momenta is closest to zero. With this solution for the x^ rest frames and 
the momentum directions qi^2 can be reconstructed, which are required for the observables 
©, ©, and dSl) . 

Fig. la shows the cos cp distributions for the production of a scalar = H and a pseudoscalar 
(|) = A in the decay channel ([5]), assuming a mass m// a = 200 GeV, both for no detector cuts 
and for applying the cuts px > 40 GeV and r\ > 2.5 (pseudorapidity) on the pions in the 
final state. In fact, the cut on r\ does not change the shape of the normalized distributions 
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shown in Fig. 1. The figure shows that the pj cuts have only a very minor influence, too. 
The slopes are given to very good approximation by the numbers below in Q. This implies 
that the shape of these distributions will be quite stable with respect to inclusion of higher 
order QCD corrections. The influence of the cuts on the shape of the distributions decreases 
for larger (|) masses. Only for light Higgs masses m^^ ^120 GeV does the chosen minimum 
Pt cut of 40 GeV have a more significant effect. 




cos{(p) Ocp 



Figure 1: Distributions of cos(p (a) and Ocp (b) for = 200 GeV and a\ polarizations LL 
or TT. 

Let us now discuss the following two situations: i) Suppose both a scalar and a pseudoscalar 
Higgs boson with (nearly) degenerate masses, for instance m// ^ ~ 200 GeV, exist and are 
produced in the reaction ([5]), / H, A. Such degenerate resonances cannot be resolved, e.g., 
in the M^x spectrum. The resulting cos cp distribution will have a shape somewhere between 
the scalar and pseudoscalar extremes shown in Fig. la, depending on the relative reaction 
rates, ii) Suppose, on the other hand, that a Higgs boson ^ with ~ 200 GeV exists^ 
which has both scalar and pseudoscalar couplings to fermions, in particular to x leptons. 
The slope of the resulting cos cp distribution will also differ from the two extremes shown in 
Fig. la. In other words, the measured distribution does not tell whether degenerate scalar 
and pseudoscalar resonances or a state of undefined CP parity were produced. This puzzle 
may be resolved using the observable Ocp- As case i) corresponds to a CP-invariant Higgs 
sector, the resulting distribution of Ocp must be symmetric (if the phase space cuts are CP- 
symmetric) and {Ocp) = 0, while case ii) will produce an asymmetric distribution and a non- 
zero average. This is shown in Fig. lb, where case ii) is illustrated with an "ideal mixture" 
(label CPmix), i.e., a (|) boson with scalar and pseudoscalar couplings of equal magnitude - 
we put ai = —bi- Again, the applied cuts have only a minor influence on the distributions. As 
already mentioned above the distributions in Figs, la, b do not change if both intermediate 

^Of course, Higgs bosons of this type might also be degenerate; for simpHcity we do not consider this 
possibility here. 
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ai are transversely polarized, while they are reflected with respect to the vertical line passing 
the abscissa value zero in the case of mixed polarizations. 

An important question is how robust is the discriminating power of these distributions with 
respect to experimental errors. In order to study this issue, using Monte Carlo methods, we 
have accounted for the expected measurement uncertainties by "smearing" the relevant quan- 
tities with a Gaussian according to exp(— (x/a)^/2), where x denotes the generated quantity 
(position in x space, momentum component, energy) and a its expected standard deviation 
(s.d.). We use here of = 15//m, of. = 15 /urn, of = 500 /urn, a^^ = O.Smrad, dE/E = 2%, 
where of denotes the s.d. of the position of the T production vertex along the beam axis, 
while of and of- are the s.d. of the positions of the respective T decay vertex along the T- 
jet axis (i.e., the direction of ai) and in the plane transverse to this axis. Furthermore, the 
uncertainty in determining the direction of ai is parameterized by an angle with s.d. Og^ 
and dE/E denotes the relative error of determining the energy of ai. These values appear 
to be realistic for the LHC experiments [19, 20]. For this simulation we use, in the case of 

= 200 GeV, a constant T flight length of 4.5 mm. 
The effect of these uncertainties on the distributions of cos(p and Ocp is shown in Figs. 2a, b, 
using again Higgs boson masses = 200 GeV. For the above set of uncertainties, scalar and 
pseudoscalar states are still clearly distinguishable (Fig. 2a), and likewise, CP-conserving 
and CP-violating states (Fig. 2b). We have made a systematic study by varying i) the masses 




-1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 -1 -0.8 -0.6 -0.4 -0.2 0.2 0.4 0.6 0.8 1 

COS(^) Ocp 

Figure 2: Distributions of coscp (a) and Ocp (b) taking into account measurement uncertain- 
ties and cuts. 



m(|) of the various types of Higgs bosons between 120 GeV and 500 GeV and ii) the expected 
measurement errors. Varying we found that the discriminating power of these distribu- 
tions does not decrease for heavy Higgs bosons. This can be understood as follows. When 
m(|) is increased the angle between the x and a[ directions in the laboratory frame decreases, 
which implies that the dependence on the smearing parameters of the distributions is be- 
coming stronger. On the other hand, the fact that the average flight length of the x leptons 
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is becoming larger reduces the sensitivity to the smearing parameters, leaving the overall 
dependence of the distributions on the above set of uncertainties rather stable. Concerning 
measurement errors we found that it is important to have under control the transverse un- 
certainty of. in the reconstruction of the x decay vertices and also the uncertainties of and 
o'q of the position of the x production vertex and of 0. In order to make use of the discrim- 
inating power of the above observables, one should achieve of. < IS/Jin, of < 30/jm, and 
Og' < 1 mrad in future experiments. Least critical are the resolution of the longitudinal x-jet 
axes (ofj and the energy uncertainty of the ai meson. Details of our results will be given 
elsewhere [21]. 

Finally, we estimate how many events ([5]) are needed in order to discriminate between i) 
a scalar and pseudoscalar Higgs boson and/or ii) between CP-conserving and CP-violating 
states, assuming = 200 GeV. As to i), we define an asymmetry A(p = [A^(coscp > 0) — 
N{cos(p < 0)]/[A^> +A^<]. From Fig. 2a we obtain from the smeared distributions = 
—0.19 and = 0.17. Thus, for distinguishing H from A with 3 s.d. significance requires 
69 events Q. Concerning ii), the result of Fig. 2b implies that for an ideal CP mixture the 
CP asymmetry Q takes the value A{Ocp) = 0.23 while it is zero for pure H, A, and degen- 
erate H and A intermediate states. Thus, 170 events ^ will be needed to establish this CPV 
effect at the 3 s.d. level. This may be feasible, depending on the masses and couplings of (|), 
after several years of high luminosity runs at the LHC [17, 18]. The x data sample will be 
increased by an order of magnitude if the above observables can be applied also to one-prong 
hadronic x decays. This is currently being investigated, employing appropriately constructed 
pseudo rest-frames [21]. 

Conclusions: The x decay channel is clearly most suited to explore the CP nature of a 
light or heavy neutral Higgs boson (|). This is an important physics issue if Higgs bosons 
are discovered. We have discussed a set of observables that serve this purpose, and we have 
shown, for Higgs boson production at the LHC and its decay via x leptons into ai mesons, 
that the above correlations and asymmetries provide powerful tools for discriminating be- 
tween CP-even and -odd Higgs bosons and for searches for CP violation in the Higgs sector. 
The measurement of these observables is challenging, but our analysis indicates that it should 
be feasible in the long run, provided enough (|) — xx events will be recorded at the LHC. 

Acknowledgements: We thank P. Sauerland and A. Stahl for helpful discussions. This 
work was supported by Deutsche Forschungsgemeinschaft SFB/TR9. 
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